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Turbulence statistics of natural airflow within a large Open Top Chamber 
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Abstract: The turbulence statistics parameters (variables) of natural airflow within a large Open Top Chamber (OTC), 4 m in high height 
and 3 m in diameter, were measured with a three-dimensional ultrasonic anemometer/ thermometer at Research Station of Changbai Moun¬ 
tain Forest Ecosystems Jilin Province, China in May 2004, for improving the field application of OTP. Results showed that because of the 
physical limitation, turbulence within OTC exhibited unique map compared with that of natural environments. There were clear daily pat¬ 
terns for most parameters. Turbulence here seemed to be isotropic and closely linked at all directions. Shape of eddies looked like a 'cylin¬ 
der' Which was very similar to the shape of OTC. Continuous airflow was always interpreted by large scale eddies from top of OTC and 
showed high interactive intermittence at all directions. 
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Introduction 

At present, there are three options for researchers to investi¬ 
gate response of plants to climate change by manipulative ex¬ 
periments, namely open-field C0 2 enrichment (FACE) system 
(McLeod et al. 1988), fully-enclosed chamber systems (Lucas et 
al. 1987) and open-top chamber systems (OTC), (Heagle et al. 
1973). Among these techniques, open-top chambers have been 
most widely used without high cost and complex environmental 
controls of closed chambers (Jones et al. 1984) or high cost of 
C0 2 required by free air release. 

Compared with FACE systems, open-top chambers alter the 
microclimate within them (Leadley et al. 1993). Measurements 
of micro environmental variables, e.g. temperatures, humidity 
and radiation, etc., with gradient ways show that there are gradi¬ 
ent changes in all factors across the chamber (Wang et al. 1994; 
De et al. 1998) which may own in large part to the altered ex¬ 
change properties, compared with the natural environment (Ham 
et al. 1991). Measurement and control of the key physical vari¬ 
ables is complex and conflict. The detailed studies on fluid dy¬ 
namics and exchange properties, particularly turbulence maps, 
will help to minimize this problem and to improve the perform¬ 
ance of open-top chamber as a field experimental tool (Leadkey 
et al. 1993). However, very few studies on these processes within 
OTC have been reported. 

On the basis of interior exchange data by using a three-di¬ 
mensional ultrasonic anemometer/thermometer in May 2004, this 
paper intends to present the detailed statistics of natural turbu¬ 
lence within a large open-top chamber and to provide basic in¬ 
formation of improving field application of OTC. 
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Methods 


Site and data collection 

Measurements were conducted in a large OTC situated at the 
Research Station of Changbai Mountain Forest Ecosystems, 
Chinese Academy of Sciences (42°24’N, 128°28’E, 738 m in 
elevation; Antu County, Jilin Province, China). The chamber, 4 
m in height and 3 m in diameter, is framed with steel structure 
and covered with glass wall. Five-year-old Pinus koraiensis 
seedlings, with an average canopy height of about 30 cm, were 
planted within the chamber. No CO z treatment was made in the 
measured chamber. A 2.5-m-tall steel pole was installed at the 
centre of OTC bottom to support the OTC. A three- dimensional 
(u, v, w) sonic anemometer (CSAT32, Campbell Inc., USA) was 
mounted at 2 m height on the steel pole, and a sensor head was 
installed at 0.25 m away from the pole. Data was continuously 
sampled at 10 Hz and recorded by a laptop computer from Au¬ 
gust First to August Twentieth, 2004. Raw data were sliced into 
30-min segments for subsequent analysis. 

Data analysis 

Linear trends were removed before the following calculations. 
Local atmospheric stability was characterized by the £ = zl L , 
where, L is Monin-Obukov length. 


L = 


6 A u' w' + w v' 


kg w'G\ 


( 1 ) 


where, k is Von Karman constant (assumed=0.4), g is the accel¬ 
eration velocity due to gravity, //is the potential temperature, and 
u, v, and w components of wind in this paper. 

By substituting u, v and wfor x in following equations, standard 
deviations of velocity components were calculated as 



( 2 ) 


Skewnesses of velocity components, Sk x were calculated as 
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Sk x = ^/a] (3) 

Kurtosis of velocity components, Kr x , were calculated as 

& t =y 4 A: (4 ) 

Integral time scales rof velocity components, were estimated as 

=r R ^) d $ (5) 

where, R x (& is the eulerian autocorrelation function at time delay 

£ 

Taylor Micro time scales. A, were estimated as coefficient of 
following regression formula, 

rM)= i-^ 2 A 2 v (6) 

Results 

Standard deviation of wind components 

Standard deviations of all wind components showed clear di¬ 
urnal patterns (Fig. 1). Daytime values were higher than night¬ 
time ones. Average values (Table 1) were lower than the reported 
values of natural environment (Baldocchi et al. 1988; Raupach et 
al. 1996; Zhang 2002). It implies that turbulence density within 
OTC is far smaller than that in natural environments. 

In contrast to the nature environments where a w is much less 
than o u and (T v (Baldocchi et al. 1988; Amiro 1990; Raupach et 
al:, 1996; Zhang 2002), standard deviations of all wind compo¬ 
nents here were very close to each other. 


Table 1. Statistics of standard deviation of wind components 


Items 

Average 

Variance 

Max. 

Min. 

O u 

0.09 

0.06 

0.26 

0.02 

Oy 

0.10 

0.06 

0.31 

0.03 

<7* 

0.10 

0.06 

0.30 

0.03 

OgOy 

0.90 

0.10 

1.15 

0.65 

OjOy, 

0.90 

0.11 

1.24 

0.56 

OytOw 

1.01 

0.12 

1.58 

0.62 


Notes: o u . a, and <r„ arc standard deviation of wind components u, v and w 



Fig. 1 Averaged daily course of standard deviation of wind components 
within OTC from 1 Aug. to 20 Aug., 2004. 

Time scales and length scales 

Taylor micro scale (A) characters contribution of turbulence 
kinetic energy (TKE) to enstrophy (Hu 1995) and size of mini¬ 
mum eddy transfer energy from the inertial subrange to the dis¬ 


sipation range. Within OTC, daytime (A u y ) was higher than 
nighttime ones (Table 2), while nighttime (X w ) was almost equal 
to its daytime value. Horizontal advection was restricted by the 
chamber and turbulence was dominant process for horizontal 
transport. Daily course of horizontal turbulence density may be 
responsible for the daily pattern of 

Integral time scales and length scales characterize the en¬ 
ergy-containing range of eddy time and spatial scales. The very 
low mean velocities within OTC imposed conditions where Tay¬ 
lor's 'Frozen eddy' hypothesis does not hold, and length scales 
based on u are no more meaningful than those based on the 
standard deviations (Amiro 1990). Following Amiro (1990), 
Length scales (L u v w ) were defined as the product of turbulence 
density (<T„. Klv ) and time scales (r„. ,,„,). 

According Table 2, at daytime, time scale and length scale of 
vertical component were much larger than that of horizontal 
components. Eddy looked like a 'cylinder' and was very similar 
to the shape of OTC. 

During nighttime, time scales of horizontal components were 
much larger than that of daytime, but for the reduction of turbu¬ 
lence density during night, the length scales increase not much as 
time scales. 

Although there was not obvious difference between nighttime 
value ( T w ) and daytime value, the nighttime value (Z^) was 
much lower than the daytime values for the big reduction of <r w . 


Table 2. Averaged Taylor micro time scales, integral time scales and 
length scales within a large OTC . 



4(s) 


r(s) 

Urn) 

Items " 

Day 

Night 

Day 

Night 

Day 

Night 

u 

0.019 

0.008 

3.08 

18.44 

0.44 

0.83 

V 

0.017 

0.008 

2.92 

19.56 

0.45 

0.85 

w 

0.005 

0.004 

17.70 

20.04 

2.98 

111 


Skewness and Kurtosis 

Skewness expresses the degree of asymmetry of a probability 
distribution. As shown in Fig. 2a, there was obvious daily pattern 
for Sk^ though the magnitude was little. Sky, for most runs were 
negative. There were not obviously daily pattern for Sk aw while 
values were almost positive. This was very similar to the airflow 
around forest canopy where u, v was positively and w was nega¬ 
tively skewed by large-scale downward movements (Shaw and 
Seginer ct al. 1987). 



Fig. 2 Averaged daily course of Skewness (a) and Kurtosis (b) of wind 
components within OTC from 1 Aug. to 20 Aug., 2004 
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Kurtosis is a measurement of peakness or flatness of a prob¬ 
ability distribution, with a value of 3 for Gaussian distributions. 
This parameter can be used to indicate the intermittent of airflow 
(Shaw 1985). Fig. 2b showed the difference between Kurtosis 
values for Gaussian distribution and that of u, v, and w velocity 
components. The Kurtosis values of all components for all runs 
were over 3. Kurtosis of wind components were highly corre¬ 
lated (Table 3). The airflow was characterized by high interactive 
intermittence at all directions. 


Equations 

A 2 

N 

P 

Kr a = 0.688 Ar,+0.644 

0.418 

648 

<0.001 

Ar u = 0.510 Kr w + 0.911 

0.237 

648 

<0.001 

Art, = 0.495 & w + 0.639 

0.253 

648 

<0.001 


Notes: Kr a , Kn, Kr w are kurotisis of wind components u,v and w. 


Summary and discussion 

We presented the turbulence statistics parameters within a 
Large OTC measured with a three-dimensional ultrasonic ane¬ 
mometer/ thermometer in May 2004. This will provide basic 
information to improve field application of OTC. 

Map of basic turbulence variables 

For the physical limitation, turbulence within OTC exhibited 
unique map, compared with those of natural environments. There 
were clear daily patterns for most parameters. Turbulence here 
seemed to be isotropic and closely linked at all directions. Spatial 
structure of eddies looked like a 'cylinder' and were very similar 
to the shape of OTC. Continuous airflow was always interpreted 
by large scale eddies from top of OTC and showed high interac¬ 
tive intermittence at all directions. 


Table 4. Statistics of Local atmospheric stability (Q within OTC 

Terms _Dav_ Night 


Mean 


- 2.02 


-6.50 



-0 20 


- 0.10 


0.00 

IV. (ms') 


0.10 


0.20 


Fig. 3 Correlation between the velocity scales for buoyancy ( W-) and 
standard deviation of vertical wind speed (a„) 


Mechanisms for turbulence within OTC 

In general, buoyancy and mechanical shear are dominant 
mechanisms for turbulence production and maintenance at local 
scales (Stull 1988). The stability indicator (Q , which characters 
the ratio of buoyancy to mechanical shear terms of TKE budget 
equation, was greater than 1 and negatively skewed all-day (Ta¬ 
ble 4). Daytime and nighttime percentages of data points for 


ICI>1 were 67 and 82, respectively. These suggest that buoyancy 
was the dominant mechanism for turbulence within OTC. Close 
correlation between the velocity scales for buoyancy ( W*), the 
velocity scales for buoyancy, and standard deviation of vertical 
wind speed (c w ) also supported this inference (Fig. 3). 

Manipulative experiments are absolutely necessarily to predict 
the behavior of terrestrial ecosystems under the changing climate. 
Open-top chambers have been most widely used without high 
cost or complex environmental controls. The basic information 
by this study and coming measurements for turbulent properties 
at multipoint within chambers will help to minimize the problem 
of measurement and control of the key physical variables in its 
field applications. 
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